Deriving the existence of BB* bound states from the X(3872) and Heavy Quark 

Symmetry 
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We discuss the possibility and the description of bound states between B and B* mesons. We 
argue that the existence of such a bound state can be deduced from (i) the weakly bound X(3872) 
state, (ii) certain assumptions about the short range dynamics of the DD* system and (iii) heavy 
quark symmetry. From these assumptions the binding energy of the possible BB* bound states is 
determined, first in a theory containing only contact interactions which serves as a straightforward 
illustration of the method, and then the effects of including the one pion exchange potential are 
discussed. In this latter case three isoscalar states are predicted: a positive and negative C-parity 
^Si — ^Di state with a binding energy of 20MeV and 6MeV below threshold respectively, and a 
positive C-parity ^Pq shallow state located almost at the BB* threshold. However, large uncertain- 
ties are generated as a consequence of the l/mg corrections from heavy quark symmetry. Finally, 
the newly discovered isovector Zi,{10610) state can be easily accommodated within the present 
framework by a minor modification of the short range dynamics. 

PACS numbers: 03.65.Ge,13.75.Lb,14.40.Lb,14.40.Nb,14.40Rt 



I. INTRODUCTION 

The existence of heavy meson bound states, a possibil- 
ity first theorized by Voloshin and Okun [1], is grounded 
on the observation that meson-exchange forces can arise 
between two heavy mesons as a consequence of their light 
quark content. In analogy with the nuclear forces in the 
two-nucleon system, these exchange forces may eventu- 
ally be strong enough as to generate bound states and 
resonances in two-meson systems. In particular, we ex- 
pect the one pion exchange (OPE) potential to drive the 
long range behaviour of two-meson systems and therefore 
to play a very special role in the formation and the de- 
scription of their corresponding bound states. The semi- 
nal works of Tornqvist 0, Q , Manohar and Wise 4] and 
Ericson and Karl Q, which already stressed the impor- 
tance of pion exchanges in the formation of deuteron-like 
meson bound states, indicated that these states are more 
probable (i.e. more bound) in the bottom sector than in 
the charm one. 

The discovery of the X(3872) state by the Belle col- 
laboration dl, later confirmed by CDF 0, DO [3 and 
BABAR probably represents the most obvious can- 
didate for a heavy meson bound state. With a mass of 
mx = 3871. 56 _± 0.22 MeV [l0| the X(3872) lies very 
close to the threshold, moo + m^.o = 3871.79 ± 

0.21 MeV, suggesting the natural interpretation of a shal- 
low bound or virtual state. The avera ge separation be- 

10 fm, will en- 



tween the and D*" mesons, 
hance strongly the role of the molecular component at 
low energies in comparison with the more compact mul- 
tiquark components. It should be noted however that the 
molecular interpretation is only possible if the quantum 
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numbers of the AT (3872) are J^"^ ^ i++ 
patible with available experimental information 
In this regard, the new experimental analysis of Ref. [13| , 
which prefers the quantum numbers 2 ^ , raises the like- 
lihood of the tetraquark and charmonium interpretations 
over the molecular hypothesis. On a related note, the re- 
cent finding of the Zfc(10610) and Zb(10650) states by the 
Belle collaboration [ij], very close to the BB* and B*B* 
thresholds respectively, might also prove to be heavy me- 
son bound states. 

The molecular descriptions of the X(3872) can be clas- 
sified in two broad types: (i) the short-range interaction 
picture [UMi, in which the X(3872) i_s the result of 
contact interactions between the D and D* mesons, and 
(ii) the potential picture [H-ii], in which the X(3872) 
is bound as the result of meson-exchanges. A middle 
path is provided by X-EFT [29*1, in which the X(3872) is 
bound due to the effect of the short range dynamics, but 
pions are included as perturbations. 

The physical motivation behind the short-range de- 
scription of the X{3872) is that the two particle conform- 
ing a low-lying bound state are too far apart from each 
other as to distinguish the specific details of the interac- 
tion binding them, a phenomenon called universality [l3l . 
However, although universality demystifies the success of 
contact theories, there are observables which cannot be 
explained without the explicitly consideration of shorter 
range components of the X(3872). A clarifying example 
is the large branching ratio |ll| . 

^^(^;i4^:^^ = i.o±o.4±o.3, (1) 

Br{X J/tpTT+TT-) ' ^ ' 

which finds a most natural explanation by considering 
the charged D+D*- component in the X(3872) [20l,|2l|- 
The potential description of the X(3872), which was 
pioneered by Tornqvist [3| , is constructed in analogy with 
the traditional meson theory of the nuclear forces [30, [HI . 
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In this picture, the D and D* mesons interact via a series 
of exchanges of increasingly heavy mesons, conforming a 
long range potential. The pion is always included, while 
more sophisticated descriptions also consider the effect of 
two pion ^28] and heavier meson exchanges (25,, 231 ■ How- 
ever, pseudoscalar and vector meson exchanges give rise 
to a singular tensor interaction, which diverges as 
at short enough distances. This singularity requires the 
regularization of the potential, usually by the inclusion 
of a form factor. Usually, these approaches do not only 
consider the DD* system, but also other heavy meson 
systems, in particular BB* . 

The conclusions which are derived from the potential 
picture are easy to summarize: in first place, meson ex- 
change forces, specially the pion, are too weak to be able 
to bind the D and D* mesons together (unless the cut-off 
in the form factor takes very large values) ; second, if the 
argument is extended to the B and B* meson system, 
the conclusion is that they are probably bound. How- 
ever, it is very difficult to determine the binding energy 
of this two-body system. This is just a restatement of 
the argument of Ericson and Karl [5] : a system with a 
higher reduced mass is more likely to bind, as a result of 
the lower kinetic energy, meaning that the existence of a 
bound state between B mesons is highly likely. 

In this work we employ two different effective field the- 
ory (EFT) descriptions of heavy meson systems at lowest 
order: (i) a pionless theory, in which the mesons inter- 
act through local contact interactions, and (ii) a pionfull 
theory, in which we include the pion exchanges explicitly 
3^S|). The EFT approach we follow 
37| , in which the potential is expanded 



(for reviews, see 
is cut-off EFT [36 

according to naive dimensional analysis 

^EFt(<7) = V^^'^q) + V^^\q) + 0(P3) ^ 



(2) 



where P represents the low energy scales of the system, 
in this case the pion mass and the relative momenta 
of the heavy mesons. Finally, a cut-off is included in 
the computation to regularize the ultraviolet divergences. 
The resulting potential is iterated in the Schrodinger or 
Lippmann-Schwinger equation in order to compute the 
observables. If pions are included, the approach basi- 
cally coincides with the Weinberg counting [s^, [s^ as 
applied in the two nucleon system, which is known to 
be very successful phenomenologically [13, |4l[ . The op- 
timal cut-off choice is to be taken of the order of the 



high energy scale of the system [42|, that is, the scale at 
which the short range effects not explicitly taken into ac- 
count in the EFT start to manifest. This corresponds, for 
example, with the energy scale at which the composite 
nature of the heavy mesons can be probed or for which 
the exchange of heavier mesons (like the p) needs to be 
taken into account, suggesting a value for the hard scale 
of Aq ~ 0.5 — 1 GeV. In this work, the natural value of 
cut-off will be roughly interpreted as the inverse size of 
the heavy mesons. 

As a consequence of the scarce phenomenological input 
available about the DD* and BB* interactions, we limit 



ourselves to a leading order (LO or order P") EFT, in 
which there is only one low energy contact operator (Co). 
The long range interaction at LO only contains the one 
pion exchange potential, resulting in a very compact and 
simple description of the two meson system. In orthodox 
implementations of the EFT concept, the counterterm 
is a running coupling constant Co = Co(A), which can 
be determined for a given cut-off A by reproducing, for 
example, the binding energy of the system (or any other 
observable). However the previous idea, although well 
suited for the DD* case, cannot be employed in the BB* 
system as there is no experimental information available, 
posing a problem for the present EFT formulation. 

The alternative we follow is the approach of Ref. [2l[ : 
we assume the counterterm Co to be saturated by short- 
range dynamics, while the cut-off is determined by re- 
producing the binding energy in the DD* system. The 
extension of the previous scheme to the bottom sector is 
straightforward once we consider heavy quark symmetry 
(HQS) [H-lill) which allows us (i) to extrapolate the sat- 
uration condition for the counterterm from the DD* to 
the BB* system and (ii) to expect the value of the cut-off 
to be similar in both cases. With these two pieces of infor- 
mation we will be able to make predictions in the bottom 
sector. This approach, which can be labelled EFT map- 
ping or saturation method, does not correspond to the 
standard formulation of the EFT concept. However, if 
the determination of the cut-off results in a natural value, 
the saturation condition will not contradict either any of 
the usual EFT tenets. The fact is that the saturation 
of the low energy constants by short-range physics is a 
sensible prospect within the EFT context, as exemplified 
in Ref. for the two-nucleon sector. The important 
point is that saturation and HQS allow us to correlate 
the charm and bottom sectors. 

The manuscript is organized as follows: in Sectllll we 
present the formalism to describe bound states between 
heavy mesons, in Sect lllll we show the prediction for a 
BB* bound state based on a contact (i.e. pionless) the- 
ory, which serves as a straightforward exposition of our 
approach. In Sect llVI we discuss the role played by the 
inclusion of the one pion exchange potential, which im- 
plies the appearance of a new BB* P-wave bound state. 
Finally, in Sect |V] we present our conclusions. 



II. FORMALISM 

In this section we will describe the two heavy me- 
son system as a non-relativistic two-body problem in 
quantum mechanics, for which the potential between the 
heavy mesons is a well-defined object. We formulate 
the bound state problem in momentum space, where the 
wave functions are obtained by solving the Lippmann- 
Schwinger equation. We also consider in detail the ex- 
plicit treatment of tensor forces, which will appear as a 
consequence of pseudoscalar meson exchange, and cou- 
pled channels. 
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A. Bound State Equation 

The Lippmann-Schwinger equation for a bound state 
reads 



\-^b) =GoiEB)V\'^B) 



(3) 



where I^'b) represents the wave functfon of the bound 
state system, V the two-body potential and Go (E) is the 
resolvent operator, which is given by 



GoiE) 



1 



E-Ho' 



(4) 



where E' < (> 0) for bound (scattering) states. In the 
expression above, Hq is the free Hamiltonian of the two 
particle system, which is defined as 



Ho\k)^'^\k). 



(5) 



where the center-of-mass motion has been removed, fi 
is the reduced mass of the system and \k) represents 
a two-particle state with relative momentum k in the 
center-of-mass system. Projecting the previous equation 
onto plane waves, we get an explicit representation of the 
bound state equation 



2n 



d^k' 



{k\V\k')^Bik'): 



(6) 



with *s(fc) = 7^ = -2^Eb and where Eb{< 

0) is the bound state energy. Finally, the solution is 
subjected to the normalization condition 



(*i3|*s) = 



(fk 

(2^ 



l*B(fc)P = l. 



(7) 



B. Central Potential 

If we assume a central potential, the previous equation 
can be projected onto the partial wave basis, defined as 

-^/4^^|fc,^m)yW(fc), (8) 

from which we can expand the bound state wave function 
as follows 



^B{k)^VA^Y.^k,lm\<^B)Yi^{k). (9) 



Ini 



For a central potential the angular momentum of the two- 
body system is conserved and the matrix elements of the 
potential in the \k,lm) basis fulfill the relation 



{k,lm\V\k',l'm') = {k\Vi\k')5ii.5,, 



(10) 



in which and 5mm' take into account the fact that 
the potential is spherically symmetric. The relationship 



between the projected potential and the potential in the 
plane wave basis is given by 
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{k\Vi\k') = — y dkdk'Yl^[k) {k\V\k')Yim{k') , [II) 

where the result does not depend on the third component 
of the angular momentum m. 

For a central potential the bound state has well-defined 
quantum numbers I and m and the partial wave expan- 
sion of ^> B{k), Eq. ([S]), simplifies to 



*s(m(A:) = V47r^'B.Kfc)i1m(fc), (12) 

where the partial wave projected wave function B,i{k) 
does not depends on m as a consequence of Eq. (jlOp . In 
such a case, the bound state equation simplifies to 



*s,/(fc) = 

2^ /■ k''^dk' 



fc2 + ^2 



27r2 



{k\Vi\k')^B.i{k'), (13) 



where the partial wave projected wave function obeys the 
normalization condition 



k'^dk 
27r2 



(14) 



Tensor Forces 



The longest range term of the interaction between a 
heavy pseudoscalar and vector meson (both of which con- 
tain a light quark) is due to one pion exchange, which 
in turn implies the existence of a tensor component in 
the potential. This tensor force will mix channels with 
different angular momenta. However, the total angular 
momentum and the parity of the two meson system will 
be conserved. Therefore, instead of projecting into the 
|fc, Zm) partial waves, we will rather consider the follow- 
ing states 

\k^ljm) = ^ \k,lmi)\l\){l\lmi\jm) , (15) 

m.1 ,ms 

where |1A) is the intrinsic spin state of the the 
pseudoscalar-vector meson system (which can be iden- 
tified with the polarization vector of the vector meson), 
and {lXlmi\jm) is a Clebsh-Gordan coefficient. In the 
previous basis, the partial wave expansion of the plane 
wave reads 

|fc,lA) = V4;?^|fc,/jm)Zji'*^(fc), (16) 

Ijm 

where the sum over I runs from |j — 1| to |j/ + 1|, and with 

Z'^^{k)^{lXlmi\jm)Yim,Ck). (17) 

Conservation of the total angular momentum implies 
that the matrix elements of the potential are diagonal 
in j and m, that is 

{kJjm\V\k'J'fm') = {k\Vw,\k')d,,,Smm' • (18) 
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In addition, as a consequence of parity conservation, the 
matrix elements between odd I and even I' , and viceversa, 
are zero. The relationship between the projected and 
unprojected potential reads 

1 

in 



A. A' 



(*), 



X (fc,lA|T/|fc',lA')Z^^ (fc'), 



(19) 



where we have taken into account that the unprojected 
potential between the pseudoscalar and vector mesons 
depends on the polarizations A and A' of the initial and 
final states. 

Owing to the properties of the tensor force, the bound 
state between a pseudoscalar and a vector meson has 
well-defined total angular momentum and parity. In such 
a case, the partial wave projection of the bound state 
wave function can be written as 

*B torn) {k)^^Yl (fc) yij,n (k) (20) 

{Op 

with '^B,ij the projected wave function, and where the 
notation {l}^ is used to denote the set of angular mo- 
menta with parity p, and with 



(21) 



It should be noted that the projected wave function \E'B,ij 
does not depend on the third component of the total 
angular momentum (m) as a consequence of Eq. (|18p. 
The partial wave projection of the bound state equation 
reads in this case as 



*B,i,(fc) 

2fi 



k''^dk' 

27r2 



{k\Vw,\k')^B,i'Ak'), (22) 



As can be seen, parity conservation implies that channels 
with I — j are uncoupled while channels with I — j ± 1 
are coupled. Finally, the normalization condition for the 
wave function reads 



E 

{0„ 



k^dk 



*B.i,(fc)P = l 



(23) 



In general we will prefer the spectroscopic notation 
"^^^^Lj to the {pjrn) notation in this work. In the spec- 
troscopic notation, S refers to the total intrinsic spin of 
the two body system, which for the particular case of a 
and D* mesons is always S" = 1 , L is the orbital angular 
momentum and J is the total angular momentum. For 
coupled angular momentum channels, we will employ a 
dash to indicate the different angular momentum com- 
ponents. For example, the p = +1, j ^ 1 state, which 
couples the S- and D- waves, will be denoted by ^Si — ^Di; 
on the contrary, the p = — 1, j = 1 state, which only con- 
tains a P-wave, will be simply the ^Pi channel in this 
notation. 



D. Coupled Channels 

In certain cases we will need to consider the existence 
of different channels in the description of heavy meson 
bound states, mainly due to isospin breaking effects. For 
example, the Al (3872) contains a neutral and a charged 
component 

|X(3872)) = ^ [ID'^D*") - \Xa) 
v2 

+ _L [\D+D*-)-\D*+D-)] \Xc), 

(24) 

where, if isospin were conserved, we would have |Xo) = 
\Xc) for the isoscalar case (/ = 0). The previous repre- 
sentation assumes that the A^(3872) is a positive C-parity 
state0. Owing to isospin breaking effects, the interaction 
can vary slightly depending on the channel and in addi- 
tion the channels can have different kinematic thresholds 
which need to be taken into account. 

The extension to the coupled-channel case is trivial, 
and only requires to consider the existence of different 
components in the wave function 



Ei*b), 



(25) 



where we use upper indices to denote the different chan- 
nels (lower indices are reserved for angular momentum). 
For the previous wave function, the bound state equation 
reads 



(26) 



/9 



where the potential V"^ is now a matrix, and Gq{E) 
takes the form 



(27) 



with the eigenvalue of the free Hamiltonian operator 
in the plane wave basis for the a channel 



Ho\k,a) ^m\k,a) = 



Aa 



\k,a). (28) 



We employ to indicate the existence of different kine- 
matical thresholds. For example, in the X(3872) the 
charged channel {D'^D*~) lies 8.06 MeV above the neu- 
tral channel {D^D*^). However, the binding energy is 
referred relative to the neutral channel. Therefore in this 
case we take Aq = and Ap = 8.06 MeV. 



^ We remind that C-parity is a good quantum number for a meson- 
antimeson system. We are also following the negative C-parity 
convention for the vector meson, in which C\D*) = —\D*), which 
is the most natural one. 
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The previous equation can be easily projected into the 
plane wave basis |A;,a) yielding 

= 

with m%{k) = {k,a\^B) and 7^ = ~2^1o,{Eb + A„), 
where we have assumed all the channels to be below 
threshold {Ea + Aq, < 0) for simplicity. The normal- 
ization condition for the wave function is given by 

E/tSi^SWP-I- (30) 

The partial wave projection of the bound state equa- 
tion can be done as in the previous cases. We will directly 
consider the general case of an interaction containing a 
tensor operator, that is 

-fcfe2E/^(Wfi^')*L,(fc'), (31) 

la ^ J 

where the wave function is subjected to the normalization 
condition 

E /5l*B,(fc)P-l. (32) 



III. CONTACT THEORY 

In this section we propose a pionless effective field the- 
ory (EFT) description of the interaction between the D 
and D* mesons conforming the X(3872). The formula- 
tion of a suitable EFT requires the existence of a separa- 
tion of scales in the physical system under consideration. 
In the particular case of a two heavy meson system and 
for momenta which are not able to resolve the finite size 
of the heavy mesons (p < 0.5 — 1 GeV), the system can 
be described in terms of the non-relativistic D and D* 
meson fields, the pion field and the local interactions be- 
tween these fields, as far as they are compatible with 
the known symmetries of the system like parity, time re- 
versal, rotational invariance and chiral symmetry. This 
description can also be applied for the B and B* system. 

The purpose of this section is to illustrate in a sim- 
ple and amenable manner the general EFT approach fol- 
lowed in the present work. For this reason, we will inten- 
tionally ignore the pion as an explicit degree of freedom 
and employ instead a contact EFT theory at leading or- 
der (LO). If pion exchanges are weak, as happens in 
the X(3872) 29], the omission of this degree of free- 
dom is not crucial and the effect will be similar to ne- 
glecting the Coulomb interaction between charged heavy 



mesons. But if pion exchanges are relevant for the de- 
scription of the system, as probably happens in the bot- 
tom sector, the decision of ignoring them will reduce the 
range of validity of the present EFT formulation to mo- 
menta below the pion mass {p < 140 MeV). The previous 
means that, if the energy of a BB* bound state is above 
Bniax ~ '^'il'^I^BB* ^ 4 McV (with rn^ the pion mass 
and ^Ibb* the reduced mass), the predictions of the con- 
tact EFT will no longer be reliable. In such a case, the 
contact theory should be regarded as a model. 

We describe the X(3872) as a bound state of a D and 
D* mesons in a positive C-parity configuration which 
interact through a momentum- and energy-independent 
contact interaction of the type 

(fc|Fc|fc') =Co, (33) 

which is regularized with a suitable regulator function. 
However, the different masses of the neutral and charged 
meson pairs, m^jo -I- mQ^.a — 3871.79 MeV and m£)+ -I- 
mjj,.- = 3879.85 MeV, plus the low lying nature of the 
X(3872), require the treatment of the neutral {D^D*^) 
and charged [D^ D*" ) components of the X bound state 
as independent and separate channels. In the neutral- 
charged basis, the contact interaction can be expressed 
as the following 2x2 matrix 

(^iycl^')-C' (} }) , (34) 

where we have assumed that the strength of the contact 
interaction is independent on whether we have a pair of 
neutral or charged mesons. 

For determining the value of the C(P^ coun- 
terterm, we follow the formulation of Ref. ^21 ], 
where a good description of the branching ra- 
tio Br{X J/iPTT+Ti-T:^)/Br{X J/V'tt+tt") was 
achieved. In agreement with the phenomenological 
model of Rcfs. [3 S^i the previous work assumes that 
the contact interaction between the D and D* mesons 
are saturated by the D meson weak decay constant /d, 
that is0 

Co^^-^--i^. (35) 

In the language of Refs. [H, the previous condition 
is equivalent to assuming that the value of the contact 
operator is saturated by the exchange of a heavy vector 
meson in the t-channel. In Ref. [21|, the counterterm 
Cq^ is regularized with a cut-off A in momentum space 
and iterated in the Lippmann-Schwinger equation. The 
value of the cut-off is determined by fixing the binding 
energy of the DD* bound state. This scheme generates a 
small isospin violation at short distances which explains 
the previous branching ratio. 



^ Notice that Refs. "S, '49*1 employ a different normalization for 
fn , which is related to ours by a factor of \/2. 
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It should be noticed that the phenomenological model 
of Refs. [H, requires the coupling of the neutral and 
charged channels to the strange one, conformed by the 
DfD*~ mesons. In such a case, the contact interaction 
can be written as the following 3x3 matrix 



\Vc\k')^C^ 



DD* 




(36) 



However, the energy of the X(3872) state is about 
210 MeV below D+D*~ the threshold. This means 
that the wave number of the strange component in the 
X(3872) is approximately 7^ = 653 MeV, representing a 
very short range contribution to the wave function, and 
approximately of the order of the natural hard scale for 
hadronic processes. Therefore we can safely ignore this 
contribution and the related strange channel in the EFT 
formulation proposed in this workQ 

The formulation of the pionless EFT for the negative 
C-parity states can be constructed by considering the low 
energy limit of the phenomenological model of Refs. [i^, 
l49j . which predicts a contact interaction identical to the 
positive C-parity case 



{k\Vc\k')c_, = C^ 



1 1 
1 1 



(37) 



with Cq^ given by Eq. (j35l) . In principle, the previ- 
ous potential implies the existence of a negative C-parity 
partner of the X(3872). However, this conclusion de- 
pends on whether we can apply the pionless EFT of 
Ref. 21] to the C ~ ~1 case without substantial modi- 
fications. In this regard, the model of Refs. [3 113 re- 
quires the negative C-parity DD* state to couple with 
(among others) the J/vt 77 and J/^P r]' channels. The first 
lies 227 MeV below the DqDq threshold and the second 
189 MeV above. This means that the C = — 1 state can 
decay into J /'^ rj with a center-of-mass momentum of 
kj/^in — 460 MeV, and has a Jj^vl component in the 
wave function with wave number ^jj^rf — 517 MeV. In 
addition, the J/ 5* 77' contribution to the wave function 
has a remarkable tendency to decay into J/^* 77 due to 
77 — 77' mixing. The aforementioned momentum scales 
are slightly below or of the order of the typical hard 
scale of hadronic processes (Aq ~ 0.5 — 1.0 GeV), and 
may require the inclusion of the J/^', the 77 and the 
Vj as explicit degrees of freedom in a pionless EFT ap- 
plicable for the negative C-parity state. If the J ^ 77 
and 77' fields are not considered, the EFT treatment of 
this channel may not be reliable unless a cut-off below 



0.5 GeV is employed. Yet there is 



^ In particular, the strange channel probability within the X(3872) 
is completely negligible. The contribution of the strange channel 
to the binding energy of the BB* state is also quite small, as we 
will see. 



no reason to believe that the structure of the LO coun- 
terterm is similar in the negative and positive C-parity 
sectors. This observation gets support from the lack of 
clear experimental evidence for the existence of a nega- 
tive C-parity DD* state in the region where the X(3872) 
is located. 

This suspicion is confirmed by the results of Ref. [lO] , 
in which the negative C-parity DD* state either (i) dis- 
appears, if the same regulator is used as in the X(3872) 
state, or (ii) moves about Ai?cm — — « 26 MeV in the com- 
plex plane, if the regulator is modified as to keep the real 
part of the energy of the negative C-parity state slightly 
below the DqDq threshold. The previous figures serve 
as a demonstration of the limits of an EFT formulation 
containing only D and D* mesons for the C = — 1 case. 
An additional limitation is provided by the fact that the 
negative C-parity partner of the X{3872) has not been 
observed experimentally, as commented in the previous 
paragraph. That is, contrary to the C = +1 case, in 
which the explanation of the J/'^uj over J/'^p branching 
ratio provides a test of the model employed for satura- 
tion, we have no analogous external check of the model 
of Refs. [H, H^l when C = — 1. In this regard, there is no 
compelling reason to trust the saturation condition given 
by Eq. |37)) . We will therefore ignore for the moment the 
negative C-parity case. 

The extension of the previous EFT formulation to the 
bottom sector is straightforward: we expect the form of 
the contact interaction to be identical in the charm and 
bottom sectors, that is 



{k\Vc\k')=C^, 



BB* 



1 1 
1 1 



(38) 



In addition, as a consequence of the tiny mass split- 
ting between the neutral {BqBo ) and charged {B^B*~) 
channels, the coupled channel structure of the contact 
potential simplifies to 



(fc|Vc|fc')-2Q 



BB* 



(39) 



The strange channel {B°Bf) lies about 180 MeV above 
the BqBq and B~^B*~ threshold. However, as a conse- 
quence of the heavier reduced mass of the BB* system, 
the related wave number of a strange component is more 
short-ranged than the corresponding one in the charm 
sector: assuming a low lying bound state, we obtain a 
wave number of 7s ~ 980 MeV. Therefore we can safely 
ignore this degree of freedom. 

The problem in predicting the existence of bound 
states in the bottom sector is how to determine the value 
of the counterterm and the cut-off. The key observation 
in this regard is that we can make sensible estimations of 
C^^* and Ab from C^^' and A^, by invoking HQS. In 
particular, we expect that (i) the saturation condition for 
the counterterm between the B and B* mesons is given 
by 



BB' 



1 

'271' 



(40) 
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in analogy with Eq. (|35|) , and (ii) the value of the coun- 
terterms, which roughly represents the inverse size of the 
heavy mesons, or equivalently, the binding energy be- 
tween the light and heavy quark, is similar in both cases, 
modulo l/niQ corrections, which leads to 



Ar = A 



X 



(41) 



where mg is the mass of the heavy quark and A^ the 
value of the cut-off in the DD* system. Provided the pre- 
vious assumptions hold, it is trivial to determine whether 
there exists a bottom counterpart of the X(3872). 



A. Description of the X(3872) State 

We describe the X(3872) state as a DD* system with 
C-parity C = -fl. Following Ref. [2l|, we include both 
the neutral and charged components, which means that 
the wave function reads 



n/2 
1 

71 



[\D+D*- 



D*+D- 



(42) 

where we are assuming and ^'-''^'^ to be S-wave. We 
will determine the wave functions of the system by solv- 
ing a two-channel bound state equation with the follow- 
ing interaction 

{k\VS^P\k')=}{^)C^^'{K)f{^), (43) 

where a,/3 = 0,C, depending on whether we are con- 
sidering the neutral or charged channel. This contact 
interaction is equivalent to Eq. p4p . In addition, we 
have regularized the contact potential of Eq. p3)) with a 
regulator f{x) fulfilling the conditions: (i) f{x) — 1 for 
a; ^ 1 and (ii) f{x) — for x 3> 1. For simplicity we 
will use a sharp cut-off regulator, J{x) — 9{1 — x), along 
this work. 

For the previous contact potential, the bound state 
equation reads 

'^''"(k) = J^fij)FiA), (44) 
+ Tx.o A 



(45) 



with F{A) a function independent of the momentum fc, 
which is given by 



i^(A) = -Co^^*(A) 



q'^dq q_ 

.2 J^K' 



27r2 " ^A^ 
X [*^«(g) + ^-^^(5)] . (46) 

We have taken = 966.6 MeV, ^ixc = 968.7 MeV 
for the reduced mass of the neutral and charged sub- 
systems. The wave numbers are defined as = 



-2^lXo [Ex + Axo) and 72.^ = -2pLXc {Ex + ^Xc)- 
Traditionally the X(3872) bound state energy is referred 
with respect to the threshold of the neutral component. 
In accordance with this prescription, we take = 
and = (w£)+ -I- m^j.-) — (m^jo -|- m^j.o) = 8.06 MeV 
for the kinematical thresholds. We take Bx = —Ex = 
0.1 -0.6 MeV. 

The bound state equation is trivial to solve, yielding 
the solutions 



.-(.)^^^/(i)^ 



■Txc 



(47) 
(48) 



where A/" is a normalization constant which can be deter- 
mined with the normalization condition 



k^dk 



[|*^«(fc)|2 + |«'^^(fc)|2] =1 



(49) 



The eigenvalue equation, which describes the running of 
the counterterm with the cut-off A, is obtained by insert- 
ing the explicit solution to the wave functions inside the 
bound state equations, Eqs. (E)) . psj) and yielding 



1 



C^*(A) 



fdq 



2tt 



r(-) 



2fi 



9' + ^Ic 



(50) 



As previously stated, we follow the approach of 
Ref. [21| in which the counterterm is saturated by /ij. 
If the binding energy has been fixed, the saturation con- 
dition determines the value of the cut-off 



C^B\A = Ax) = -^, 



(51) 



where we take the value /d — 210 ± 10 MeV, a value 
compatible with recent lattice simulations. The satura- 
tion condition can only be fulfilled for a specific value of 
the cut-off A = Ax- For the range of bound state ener- 
gies considered in this work, and a sharp cut-off function, 
we obtain 



A 



X 



555i44MeV, 



(52) 



where the error accounts for the binding energy window 
and the uncertainty in fjj B The cut-off can be physically 
interpreted as the scale at which the system starts to 
resolve the composite nature of the D and D* mesons. 
That is, the cut-off is related to the size of the heavy 
mesons. 



'1 In contrast, Ref. Hil takes /j, = ^2 X 165 MeV (i.e. 233 MeV) 
and Bx = 0.1 MeV, for which the cut-off Ax = 653 MeV is 
obtained. 
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B. The BB* Bound State 

The description of the theoretical BB* bound state 
involves a neutral and a charged component with positive 
C-parity, C = +1, that is 

v2 

+ J- [\B+B*-) -\B*+B-)] 

(53) 

As we will see later, the isospin breaking effects are neg- 
ligible in this case, and the neutral and charged wave 
functions are approximately equal. In analogy with the 
DD* case, the BB* interaction is taken to be 



where a, /3 



^)C„^^-(A)/(^ 
0, C. The bound state equation reads 



{k\v^'^\k') = fi^)cr'wfi-) 



2n 



_BB'(0) 



^ 'SB*(0) 



(54) 

(55) 
(56) 



with given by 
F(A)=- 



(57) 



We take f^BB'fi — 2651.1 MeV and f^BB'.c = 

The wave num- 
+ Abb*,o) and 
). In this case, 



2651.01 MeV for the reduced masses 
bers are 7|s.(o) -2/^sb*(o) 

„2 



^BB'(C) 



^BB'{0) ' 



SB* 



(C) 



SB* 



^BB*(C) 



the kinematical threshold gaps are given by A^^*(q) = 
and A 



BB'{C) 



(0) 

(tob+ + mB*-) — {mBO + niB"') 



-0.35 MeV (with an error of 0.27 MeV). 

The solution of the bound state equation is analogous 
to that found in the DD* system. Taking into account 
the size of the threshold gap, A^^*^^) ~ 0.35 MeV, we 
can make the approximation 



*^^*(o)(fc) = *^^*(^)(fc). 

'^^*(fc) 



V2 



1 + 



A 



SB* 



E 



BB' 



,(58) 



where the relative error will be fairly small for bound 
state energies above 3 — 4 MeV. The isospin symmetric 
wave function is given by 



SB* 



BB' 



.1. «x) 



(59) 



where A/" is a normalization constant which may be ob- 
tained by 

Pdk 



27r2 



hi' 



SB* 



(60) 



In the isospin symmetric limit, the the eigenvalue equa- 
tion reads 



1 



'fd.'l f2f1^ '^t^BB' 



27r2 ' 'A' q^+-il ' 



yiv, J H • IBB 

and finally, the saturation condition is 



1 



C^" (A = Ab)--^, 



(61) 



(62) 



with Jb = 195 ± 10 Me_V. From HQS we expect the value 
of the cut-off in the BB* system to be similar to the DD* 
system, that is A^ ~ Ax, as already mentioned. This 
assumption is equivalent to presuming the size of the D 
and B mesons not to depend strongly on the mass of the 
heavy quark. 

In this case the known and unknown parameters are 
different than in the DD* system. Contrary to the 
X(3872) case, the bound state energy of the theoretical 
BB* state is unknown and needs to be determined. In 
this regard, we will use the eigenvalue equation Eq. (|6ip , 
the saturation condition Ec^. (15^ and the similar cut-off 
assumption to predict the expected bound state energy 
of the BB* system. 

The calculation has three error sources: (i) the bound 
state energy of the X(3872), which ranges from Bx = 
0.1 — 0.6 MeV, (ii) the uncertainties in the value of /s, 
and (iii) the error in the relation A^ ~ A^. The third 
source of error is unknown, but can be estimated from 
heavy quark symmetry: we expect the cut-offs to be 
equal modulo 1 /mQ corrections 



Ar = A; 



\mQ 



(63) 



being mg the heavy quark mass. The problem is how 
to evaluate the size of the l/mg corrections. A possible 
estimation comes from the observation that y^mjjfu and 
^wlbIb also differ by l/mq corrections 



ZmBfB = \/mDfD + O 



mq) 



(64) 



Assuming that the relative size of these corrections is the 
same for the constants and the cut-off, we obtain 



As - Ax 






As 




y/rriB /b 



- 1 



(65) 



which yields a ^ 10 — 15% relative error to the naive re- 
lation Ab = Ax- Therefore, we take A^ = Ax (1 ± 0.15) 
as a conservative estimate. For the X(3872), we ob- 
tained from the saturation condition Ax = 555^44 MeV, 
in which the uncertainties coming from Bx and fu are 
already taken into account. Considering the l/mg cor- 
rections and quadratic error propagation, the cut-off win- 
dow for the BB* raises to Ab = 555^94 MeV, where the 
errors are dominated by the l/mg term. 

For the previous estimations of /b and A^, we ob- 
tain a bound state energy of B = 45^35 MeV. The 
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binding energy is clearly beyond the limits of a contact 
range theory, which means that the inclusion of pion ex- 
changes is mandatory from the EFT viewpoint. The 
explicit inclusion of isospin breaking shifts the bound 
state energy about AS — — 0.2MeV for the central value 
B = A5 MeV, much below the error in the determination 
of the bound state energy, and is therefore a negligible 
effect. The effect of including the strange component in 
the charm and bottom sectors is small, generating a shift 
of AB = —4 MeV of the BB*, again much smaller than 
other uncertainty sources; the probability of the strange 
state in such a case is Ps — 2 — 3%. 



a I 



P 



p 



FIG. 1. Diagram representing the pion-vector-pseudoscalar- 
meson vertex. The thin line represents a pseudoscalar meson, 
the thick line a vector meson with polarization i, and the 
dashed line an incoming pion with isospin a and momentum 
9- 



IV. THE ROLE OF ONE PION EXCHANGE 

In this section we consider the role of one pion ex- 
change in the DD* and BB* systems. As in the two 
nucleon system, we expect the long range interaction be- 
tween two heavy mesons to be driven by one pion ex- 
change. From the EFT perspective, the inclusion of mul- 
tiple pion exchanges, which are in turn constrained by 
the requirement of broken chiral symmetry, extends the 
range of validity of the effective description up to mo- 
menta of the order of the p mass or the binding energy 
of the heavy and light quarks conforming the heavy me- 
son. Both scales are of the same order of magnitude, 
yielding a breakdown scale of Ag ^ 0.5 — 1 GeV. In this 
respect, we expect the EFT formulation to be able to 
describe bound states up to a maximum binding energy 
of Binax ~ Al/2^iBg, 50 - 200 MeV for the particular 
case of the BB* system. 

As already noticed in Ref . Q , an interesting difference 
with respect to the two nucleon case lies in the role of the 
short range interaction, which does not involve a strongly 
repulsive core in the heavy meson case. In fact, the short 
range force is attractive in S-waves, as can be deduced 
from the existence of the A" (3872) state. This factor will 
help the formation of S-wave bound states in the BB* 
system. 



A. The One Pion Exchange Potential 

1. Denvation of the One Pton Exchange Potential 

The one pion exchange potential between a pseu- 
doscalar (P) and a vector (P*) meson (both containing 
a light quark) can be derived from the PP*7r vertex of 
Fig. ([T]). In Heavy Hadron Chiral Effective Field The- 
ory [50, 51], the non-relativistic amplitude corresponding 
to this vertex in the isospin basis is given by 

A[P*'^^PT,-)^j-^e,-q (66) 

with Ta the isospin operator (i.e. the Pauli matrices) for 
the pion 7r° in the cartesian basis and ti the polariza- 



tion vector of the P*^ meson. For the pion weak decay 
constant, we use the normalization /^r ~ 130 MeV. 

The value of g represents the coupling of the pion to 
the light quark in the pseudoscalar/ vector meson. In non- 
relativistic quark models, we expect g = 1, while in chiral 
quark models, g = 0.75 (see, for example, Ref. fs^). 
However, the best way to determine the value of g is 
from the D* decay width. For the particular cases of the 



decays D*^ D°n+ 
amplitude reduces to 

A{D*+ ^ 
A{D*+ ^ 



and D* 



ZJ+TT^, the previous 



D 



V2U 



(67) 
(68) 



as can be trivially checked, yielding the well-known decay 
rates 



T{D* 



T{D* 



9 



= 9 



(69) 



(70) 



From these decays (and taking into account the addi- 
tional electromagnetic decay D*"*" — > I?+ 7), Refs. [sslls^ 
obtain the value g = 0.59 ± 0.01 ± 0.07. For simplicity, 
we will take g = 0.6 ± 0.1 in the DD* sector. 

For the case of the BB*Tr vertex, the value of g cannot 
be determined from pion decay. However, from heavy 
quark symmetry we should expect a value similar to that 
in the DD*Tr vertex. This expectation seems to be con- 
firmed by lattice simulations [55n58|, which suggest a 
value in the range 0.5 — 0.6. We will take g = 0.55 ± 0.10 
in the BB* sector. However there are also other deter- 
minations suggesting a smaller value of g. For example, 
Ref. m obtains g = 0.37l^;^2 

based on a study of the 
Dyson-Schwinger equations. Therefore, we will also dis- 
cuss the consequences of having smaller values of g in the 
BB* system at the end of this section. 

In the static limit, the one pion exchange potential 
in momentum space between a pseudoscalar and vector 
meson (see Fig. ^) takes the form 
■,2 



VbpE(g) 



9 



Tl ■ Tl 



(71) 



q -I- yU 



10 



I 



iViq)= J- 



q 



I 



FIG. 2. Diagram corresponding to the one pion exchange 
potential between a pseudoscalar and a vector meson. 



in the potential. This problem is overcome by taking a 
principal value prescription 



1 



= V 



1 



2q 



6iq~tl), (76) 



where p,^ — and V denotes the principal value. This 
is equivalent to the change — > i/i in the coordinate space 
potentials. We will ignore the imaginary piece of the 
potential, which is related to the decay of the X{3872) 
to DDtt. a more complete treatment of the < case, 
in which the DDtt channel is included explicitly, can be 
found in Ref. I60ll. 



2. C-parity 



where t^2) is the isospin operator in vertex 1(2), ^i[2) is 



the P* meson polarization vector and /i 



to; 



(TOp. 



nip) , with mp* and rnp the masses of the vector and 
pseudoscalar mesons. For the pion propagator we have 
used Dabiq) — Sab/iq^ — toJ — ie), where a, b are isospin 
indices. In the isospin symmetric limit, we will take for 
the pion mass the value = 138.03 MeV (which corre- 
sponds to (tott" + 2to^±)/3). The small imaginary piece 
(ie) is added for dealing with the /i^ < case. 

The form of the OPE potential in configuration space 
can be obtained by Fourier transforming the momentum 
space representation of Eq. ((7T|) . in which case we obtain 

VoPKir)^ I 7^VbpE(<7)e-'«-^^ 



(27r)3 

= ^r,.n(?/.V)(rr.V)^. (72) 
The previous expression can be rewritten as 



VbpE(f') = ■ Ti • ei 6{f 
+ T2 • fi 6*2* • eiWc{r) 



(3 6*2* • f £*! 



ei)WT{r) , (73) 



where the central and tensor components of the potential, 
Wc and Wt, are defined as follows 



2, ,3 



Weir) 
Wrir) 



~fir 



g M 

247r/2 fir 



247r/2 fir 



1 



fir /i^r^ 



(74) 
(75) 



It should be noted that by making the replacements — > 
<?ii 9 ~^ 9Ai— 1.26) and fi ttItt we recover the one 
pion exchange potential for the two-nucleon system. This 
representation will be useful for determining the channels 
in which one pion exchange is more attractive. 

A particular problem which arises with the treatment 
of the one pion exchange potential in the DD* and DD* 
systems is that fi^ is negative, generating a singularity 



We are interested in meson- antimeson systems where 
C-parity plays an important role. In this regard, it should 
first be noted that the previous OPE potential has been 
obtained for a meson-meson system. The corresponding 
potential for the meson-antimeson system in the isospin 
basis can be obtained by means of a G-parity transfor- 
mation, in which case the potential remains unchanged 

Vo''p''E(^1 = ^OPE(r), (77) 

where VbpE('^) is the OPE potential as given in Eq. (|73p . 

Now we take into account that the previous meson- 
antimeson potential is written in the following basis 

V^^^^if)^{PP*\V\P*P) 

^ {PP*\V\P*P) , (78) 

where C-parity plays no role. Therefore, we project the 
potential into a basis with well-defined C-parity, that is 

\PP*{rf))^^[\Pn-rf\P*P)] , (79) 

where r] represents the intrinsic C-parity of the system 0. 
In this basis, the potential reads 

'yin + Vout TrPp" 



{PP*{rfout)\V\PP*{Tfin)) 



(80) 

that is, the OPE potential conserves C-parity and the 
sign of the potential depends on the intrinsic C-parity of 
the system. 



Of course this depends on the C-parity convention for the heavy 
vector meson. In the present work we have taken C\P*) = 
— \P*), which also implies an extra minus sign when perform- 
ing the G-parity transformation on the vector meson. If we take 
into account the negative G-parity of the pion, the potential 
in the isospin basis does not change sign for antiparticles. If 
we had chosen the plus sign convention in the C- and G-parity 
transformations of the heavy vector mesons, the sign of the po- 
tential would have changed. However, in such a case we would 
also have obtained an additional minus sign after projecting into 
well-defined C-parity states in Eq. I I80I I. thus leaving the final 
form of the potential in the C-parity basis unchanged. 
" The total C-parity of the heavy meson-antimeson system will be 
given by C = (—I)'?}, where I is the orbital angular momentum. 
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3. Isospin Breaking Effects 

In the previous sections we have derived the potential 
in the isospin symmetric basis. However, isospin may be 
broken due to two effects: (i) the different masses of the 
charged and neutral pions and (ii) the different masses of 
the charged and neutral heavy mesons. For the X{3872) 
(the DD* system), the most important effect is (ii), due 
to the weakly bound nature of the X{3872). On the 
contrary, for the BB* system we expect both types of 
isospin breaking effects to play a relatively minor role. 

To take into account isospin breaking, we consider 
the charged and neutral components of the system sepa- 
rately, i.e. 

\PP*iv)) = «o \{PP*{v))o) + «± \{PP*{V))±) ■ (81) 
In this neutral-charged basis, we make the substitution 



(82) 



and the OPE potential reads 



''ope 



-V 



-2t4±((f,^±) ~V^o{q,fiQ) 



with {q, fJ-) given by 



F^((f,/i) = - — 



£2 ■ gei ■ g 



(83) 



(84) 



where fi must be evaluated for the particular channel un- 
der consideration. For the neutral-neutral channel we 
have Hq = m^o — {irip.a — mpo)'^, while in the charged- 



charged channel /j,'q 



{mp-± — mp±) 



For 

the term of the potential connecting the charged-neutral 
channels there are also two possible definitions for the 
value of /i-t, which are /ij. = "^^± ~ {mp,± — mpo)'^ and 
fj,'^^ = rri^j^ — {mpto — mp±)'^. They are different as a 
consequence of mp*± — mpo ^ mp,o — mp±. To include 
this correction, we simply average over the two possible 
definitions of fi± 

1 



v^± {q, M±) = ^ {K± (q, ^J■±) + K± {q: m±)) 



(85) 



4- The Relative Strength of One Pion Exchange 

The strength of the OPE potential in the PP* system 
is different in each partial wave (i.e. ^^+^Lj channel). If 



We need to take into account that the D" and D*" mesons have 
an extra minus sign when expressed as isospinors, that is and 



I 2 2'- 



we ignore the effect of zero-range contributions, bound 
states between two heavy mesons will be more probable 
in channels for which the long range potential (in this 
case, OPE) is most attractive. For this purpose we em- 
ploy the coordinate space formulation, which simplifies 
the identification of the attractive channels. In addition, 
to avoid the problems derived from the existence of a neu- 
tral and charged component of the wave function, which 
will obscure the analysis, we take the isospin symmetric 
limit in the following discussion. The methods employed 
in this subsection were originally developed in Refs. |6ll - 
|63| within the context of the renormalizability of nuclear 
forces in chiral EFT. 

The configuration space wave function of the heavy 
meson pair (with total angular momentum j) can be de- 
composed into a radial and angular piece 



B{pjm) \ 



(86) 



where uij are the reduced wave functions of the PP* 
system and {l}p, j, m and yijm were already defined in 
Sect. Ill Cl Ignoring zero range contributions, the reduced 
Schrodinger equation for the wave functions uij (r) reads 

- u'i'j{r) + 2^pp. Vwj{r) ui,j{r) 
I' 

+ = -7pp. uij{r) , (87) 

where /ipp. is the reduced mass of the PP* system, 
7pp, = fippt Eppt the wave number of the bound 
state, and Epp, the center-of-mass energy of the two me- 
son system. The partial wave projection of the potential 
can be obtained from the following expression 



E 

A, A' 



x(lA|y(rl|lA')Z^ if), 



(88) 



as can be checked by adapting the methods of Sect. Ill CI 
to coordinate space. 

For the particular case of the OPE potential, Vwj reads 

Vwjir) =-VT [Sii'Wcir) + Sw^Wrir)] , (89) 

where 77 is the intrinsic C-parity and t = T2 • Ti = 
21 {I -I- 1) — 3, with / the total isospin of the two meson 
system. Wc and Wt are the central and tensor pieces of 
the potential, see Eqs. and The Sw and Swj 

represent the partial wave projection of the 6*2* ' £1 ^^^d 
'S'i2('^) {— 3 £2* ■ rei ■ f — £2* • Cl) operators respectively. 
The S{f) contribution to the OPE potential, see Eq. (|73p . 
has been ignored as it is a zero range contribution. 

The matrix elements of the tensor operator in the par- 
tial wave basis are the following: (i) for uncoupled chan- 
nels {I = j) we have Sjjj = — 1, (ii) for coupled channels 
{I — i ±1) we obtain 



1 



2j + 1 V-3V7(7+T) 



-3^7(7+1) 



(90) 
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where the matrix notation stands for {Sj)ii' = Swj 0; 
finally, (iii) for the '^Pq channel, which is uncoupled (but 
for which I = j + 1), we have 5iio = 2. 

Naively we expect the strength of the OPE potential to 
depend on the interplay between the central and tensor 
components. However, at distances below the pion wave- 
length (fir < 1), the central piece contribution is small 
in comparison with the tensor component 0. In particu- 
lar, we expect the behaviour of the tensor force to 
determine the relative strength of OPE effects. At short 
enough distances, /ir ^ 1, we have 

2pipp, Wrir) ^ ^ , (91) 

where Rt is a length scale related to the strength of 
tensor OPE. For the OPE potential, we expect in turn 

2fipp,Vu.,ir)^ Xu'j^, (92) 

where Xwj = —rjTSwj. That is, the Swj matrix deter- 
mines which channels are attractive. 

For the uncoupled channels the identification of the 
most attractive channels is trivial as it only requires the 
multiplication of the matrix element of the tensor force 
by C and r. For the coupled channels, however, the 
relevant observation is that at short enough distances, the 
tensor force (~ l/''^) overcomes the centrifugal barrier 
(^ l/r-^) and the central force (~ 1/?"), in which case we 
can approximate the full Schrodinger equation by 

I' 

= -7pp.«y(0+o(^), (93) 

where we have ignored the centrifugal barrier and poten- 
tial contributions which diverge less strongly than 
If we define the following vector, which contains the two 
components of the wave function 

^ (5;;::) • 

we can rewrite the Schrodinger equation above as 

- <M -V-^ S,", = -7|.p. u, + O(^) . (95) 

In this notation, it is clear that any transformation which 
diagonalizes S^- also diagonalizes the Schodinger equation 
at short distances. This means that the attractive or 



As a curious fact, the matrix elements of the tensor operator for 
the PP* system have opposite sign and half the strength of the 
corresponding matrix elements in the two-nucleon system. 
Indeed this is the case in the two nucleon system, where only the 
tensor component of the OPE potential needs to be treated non- 
perturbatively, while the central piece can always be regarded as 
a perturbation. 
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+1 
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TABLE I. Relative strength of the tensor component of OPE 
for the different partial waves in the isospin symmetric basis. 
The I{rj) notation indicates the isospin and the intrinsic C- 
parity subchannel under consideration. Negative values in the 
table denote channels in which the tensor force is attractive. 
In the coupled channel case, we show the two eigenvalues of 
the tensor force. The strength of tensor OPE in the periph- 
eral uncoupled (coupled) waves coincides with that of the ^Pi 
(^51 - ^Di) channel. 



repulsive character of the tensor force is determined by 
the eigenvalues of Sj, that is 

R, S, Rj = (^"J ^) , (96) 

where Rj is the rotation matrix which brings Sj into 
the diagonal basis. Independently of the values of rj and 
T, there is always an attractive and a repulsive eigen- 
channel. In this regard, at short enough distances the 
coupled wave function contains an attractive component 
which may generate a bound state. 

The relative strength of the tensor force for the dif- 
ferent partial waves, that is, the eigenvalues of Xwj, can 
be found in Table [D The most attractive combination, 
A = —6, is achieved in isoscalar (/ — 0) channels with 
negative C-parity (77 — —1), and in particular the uncou- 
pled ^Po channel and all the coupled channels (^5*1 — ■^£'1, 
^P2 — ^F2, etc.). Next in attractiveness, the isoscalar 
(/ = 0) channels with positive C-parity (?/ — +1) are 
to be found, such as the uncoupled ^Pi and the cou- 
pled ^Si_ — ^Di. However, owing to the repulsive role of 
the centrifugal barrier, we only expect the lower partial 
waves, such as the ^Si — '^-Di, the ^Pq and the "^Pi, to be 
the most promising candidates for a bound state. 

In Table [Til we list the minimum value of the momen- 
tum space cut-off for having a bound state in the lowest 
partial waves. We consider the full OPE potential, as 
defined in Eqs. ([7T|) or ([75)) . That is, we include the 
zero-range 5 contribution. As we can see, the previous 
expectations are approximately fulfilled expect for the S- 
waves: the deviations in the ^Si — ^Di channel are caused 
by the zero range piece of the OPE potential. In the 
^P2 — channel we do not find bound states for values 
of the cut-off below 2 GeV, probably due to the strong 
centrifugal barrier of the F-wave subchannel, ruling out 
the possibility of a molecular interpretation in case the 
J^^ = 2 ^ assignment [l^l turns out to be confirmed 
by future works. 
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TABLE II. Cut-off (in MeV) at which the first bound state 
appears for the OPE potential in different partial waves. The 
I{ri) column specifies the isospin and the intrinsic C-parity 
subchannel (for the DD* , where isospin breaking has been 
taken into account, / is to be interpreted as "mostly a 7 = 0/1 
state with some small admixture of 7 = 1/0"). Partial waves 
without bound states below A = 2 GeV are not displayed. 
Errors take into account the uncertainty in the PP*tt coupling 

g- 

5. The Negative C-panty States and the P*P* System 

An interesting feature of the OPE potential, which we 
have ignored until now, is that it can mix the PP* and 
the P*P* heavy meson systems. The reason is the fol- 
lowing P*P*Tr vertex 

A{P*' -> F*^'^) = * (el X e*) ■ q, (97) 

which allows the PP* — > P*P* transition to happen via 
the one pion exchange mechanism depicted in Fig. ([3|). 
However, owing to the energy separation between the 
PP* and the P*P* thresholds, the coupling between 
these two systems will not become evident unless the en- 
ergy of the PP* bound state is comparable to the mass 
difference between the P and P* heavy mesons. In the 
bottom sector, this requires a binding energy of about 
50 MeV: for a BB* bound state at threshold, the wave 
number of the B*B* component is ^b'B' — 0.5 GeV, 
which is of the order of the hard scale of the system 
(Aq ~ 0.5 — 1.0 GeV). Consequently we can ignore the 
eventual B*B* short range component of the wave func- 
tion of a BB* bound state at the price of setting the 
hard scale at the lower bound Aq ^ 0.5 GeV E3. From 
the power counting perspective this prescription can be 
taken into account by noticing that the mixing between 




q 



FIG. 3. One pion exchange diagram mixing the pseudoscalar- 
vector (PP*) and the vector- vector {P*P*) heavy meson sys- 
tems. 



the PP* and P*P* channels is suppressed as the ratio 
between the P*P* and PP* propagators near the PP* 
threshold, that is, by two orders in the EFT expansion 

Gf^'{E) _ q^-2fipp.E fP_y 

GPP\E) q^ + 2^ipp,A-2^ipp,E^\Ao) ' 

(98) 

where q ^ ^J2[ipp, E ^ P are considered to be light 
scales, while yj2^pp, A ~ Aq is a heavy scale. In this 
regard, we expect the mixing of the channels to be an 
effect similar in size to chiral two pion exchange (TPE), 
which also enters at 0{P'^). 

If we are interested in a more detailed account of the 
effect of the P*P* channel in the description of the PP* 
bound states, it should be noted that the conservation 
of parity, C-parity and total angular momentum requires 
that the coupling between the PP* and P*P* channels 
only occurs between states with the same J^'^ quantum 
numbers. In particular the intrinsic negative C-parity 
PP* states will tend to have a stronger mixing with the 
P*P* system than the rj = +1 states. The reason is that 
the '^^'^^ Lj{PP*) partial wave with rj = —1 has the same 
quantum numbers as the '^^'^^Lj{P* P*) wave!"]. In this 
regard, the J^*^ = 1"' PP* state is expected to have the 
strongest mixing with the P*P* system, as the coupling 
can happen via S-wave_(3S'i(PP*) - ^Si_{P*P*)). In the 
0"+ states, the ^PoiPP*) and ^Po{P*P*) channels can 
also mix, although being a P-wave the ^Po{P* P*) com- 
ponent is expected to be suppressed. The overall effect 
of this mixing will be attractive: the short range P*P* 
component of the wave function can always take a config- 
uration which minimizes the energy of the system, thus 
increasing the binding. 




The situation is in fact very similar to ignoring the A isobar res- 
onance in the low energy description of the two nucleon system, 
where a similar wave number is obtained for the eventual A^A 
short range component of the deuteron. 



We remind that the C-parity of a PP* and p* p* systems is 
C = {—IYt] and C = (—1)'+" respectively. This implies that 
intrinsic negative C-parity ^^'^^ Lj{PP*) states can always mix 
with the 2'S+l£,j(p*p*) states. 
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On the other hand, the intrinsic positive C-parity PP* 
states will tend to couple much more weakly with the 
P*P* channel. The extreme case is the ^Po{PP*) partial 
wave {J^'~^ = ), which does not contain any P*P* 
component as there is no P* p* state with such quantum 
numbers. In the 1++ state, the mixing is expected to 
be small: the ^Si — ^Di{PP*) channel can only couple 
with the ^ Di{P* P*) partial wave. In these cases we can 
probably ignore the PP* mixing without significantly re- 
ducing the convergence of the theory. 

B. The X(3872) with One Pion Exchange 



where a, (3 = 0,C and = Vq^e + • We will 

employ a sharp cut-off regulator, f{x) = 9{1 — x). 

As in the contact theory case, the value of the cut-off is 
determined from the saturation condition, Eq. (|5ip. For 
the central values g = 0.6, Bx = 0.35 MeV and fn — 
210 MeV, and a sharp cut-off regulator, the saturation 
cut-off 

Ax = 396^^1 MeV (104) 

is obtained, where the spread in the cut-off value comes 
from the uncertainty in g, Bx and Jd E3- 



We describe the X{3872) as a ^^i-^Di, C = +1 state, 
where the charged and neutral components are treated 
independently. The wave function reads 



+ 



1 
1 

72 



7=)*0\ 



ID'-'D 



D*^D°) \ ^ 



(pjm) 



{k) 



\D+D*-) - \D*+D- 



(jpjrr 



(99) 



where p = +1, j = 1 and {l}^ = 0,2. The *fp°„,)(A:) 

and ^^^„j') {k) wave functions can be decomposed into 
an S-wave and D-wave component 



(+1 Im) 



(fc) 



*s°(fc)yoi™(fc) 



+ *:^°(fc)3^2i™(fc) , (100) 



where the subscript a = 0, C denotes the neutral and 
charged components, and where we have used S and D 
instead of Z = 0, 2 for concreteness. The normalization 
condition for the wave function reads 



E 

Xa. 



k^dk 



l*s°(fc)r + l*D (fc)r =1- (101) 



The wave functions are obtained by solving numer- 
ically a four channel Lippmann-Schwinger equation, 
Eq. (PT|) . with the OPE potential and a contact term. 
The OPE potential in the neutral-charged basis is taken 
from Eqs. ((831) and ([84]). For < 0, we employ the prin- 
cipal value prescription and ignore the imaginary piece 
of the OPE potential. The partial wave projection of 
the OPE potential can be obtained from Eq. ^T9\\ . For 
the contact piece, the partial wave projection is trivial, 
yielding 



{pJj\vf\p',l'j)^C^'''Si,oSv,o 



DD* 



(102) 



The full potential is regularized with a suitable regulator 
function f{x) 



(pJjlV^f^lp'J'j) ^ f{'^){k,lj\V'^^\p',l'j)f{^). 



(103) 



1. The Negative C-Parity State 

An interesting advantage of including pion exchanges 
explicitly is that the EFT framework can be now applied 
to the negative C-parity sector. With OPE the satura- 
tion cut-off takes a smaller size than the typical momenta 
associate to the J/'^rj and J/'^rj' channels, which means 
that we are within the range of applicability of the the- 
ory. In this regard, while the contact theory seems to 
predict a negative C-parity partner of the X(3872), the 
inclusion of OPE in the I{J^^) = 0{V ) channel pre- 
vents the formation of a bound state unless the cut-off 
A is greater than 760 MeV. However Ax turns out to 
be too small, meaning that OPE explains in a natural 
manner why the negative C-parity state does not exist. 



C. BB* Bound States with One Pion Exchange 

The number of bound states between the BB* mesons 
depends on the natural value of the cut-off for this sys- 
tem. Assuming the relation A^ ~ Ax, and taking into 
account the results from Table ITIl we expect to find two 
states in the BB* system: a I{J^'^) = 0(1++) and a 
/( J^C) ^ 0(0-+) state. The 0(1++), - ^Di statejs 
almost universally predicted in works investigating BB* 
molecular states, see for example Refs. 19, 25-27]. The 
prediction of a ^Pq resonant state is less usual, as most 
works concentrate in the S-wave states. In addition, ow- 
ing to the fundamental role of contact interactions in 
the formation of heavy meson bound states, we will find 
that the negative C-parity 0(1^ ) state is bound too. In 
Table IIIII there is a summary of the binding energies of 
the bound states and resonances predicted in the present 
work. 



In particular, we have Ax 



396 



+38 +7 +19 
-37 -10 -19 



MeV, where the 



errors come from g, Bx and respectively. For the final value, 
we have assumed that the errors add quadratically. 
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OPE 
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fl 7+00 1 -1 +00 i 
'J-'-l.D ^-^-1.62 
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TABLE III. Summary of bound states and resonances in the 
bottom sector. For the bound states we show the binding 
energy {B = —Ecm), while for the '^Pq resonance we indicate 
the center of mass energy of the state. 



1. The Positive C-parity ^Si — ^Di Bound State 

In first place, we consider a theoretical BB* bound 
state in the '^Si — '^Di channel with positive C-parity, 
total isospin / = and J^^ = 1++. For this system, the 
wave function reads 



1 



(pjm) (^') - 2 



\B"B*") - \B*°B°) 



\B+B*-) - \B*+B-) 



(pjm) 

(105) 



where p = +1, j — 1, {l}^ = 0,2 and we have assumed 
isospin symmetry to hold. As in the X(3872) case, the 
wave function can be decomposed into an S- and D-wave 
piece 



(106) 



*( + llm)(fc) = V47r ^ 4's(fc)3^oim(fc) 

+ *I3(fc)J^21™(fc) 

subjected to the normalization condition 
k^dk 



27r2 



[\^s{kr + \^D{k)\' 



1. 



(107) 



For the BB* system, as a consequence of isospin sym- 
metry, we solve a two channel Lippmann-Schwinger equa- 
tion, Eq. pip . The full potential results from adding the 
OPE contribution, as described in Eq. d7T]) , and the con- 
tact theory contribution, which in the partial wave basis 
reads 



{p, lj\Vc\p',l'j) 



(108) 



where a factor of 2 needs to be included as we are working 
in the isospin symmetric basis. We assume the contact 
interaction to be determined by the saturation condition, 
Eq. dSS]), where we take /s = 195 ± lOMeV. The fuh, 
partial wave projected potential is then regulated with a 
regulator function f{x), 



{pJj\V\p'J'j) 



f{j){k,lj\V\p\l'j)f{'j^] 



with V = Vqpb + Vc- As in previous cases, we will em- 
ploy a sharp cut-off, f{x) — 9{1 — x). The cut-off is sup- 
posed to be As = Ax, modulo l/mg corrections, which 
we estimate to be of the order of 15%. This translates 
into the value A^ = SOGl^gMeV, where the combined 
uncertainty is the result of quadratic error propagation. 
The cut-off error is dominated by the estimation of the 
size of the l/mg corrections. 

For the values of g, /s and A^ we are using, we ob- 
tain an estimation of the bound state energy oi B = 
20l]^2 MeV. The error estimation is completely domi- 
nated by the cut-off uncertainty, which alone would gen- 
erate B — 20't_\i MeV. The other sources of error estima- 
tion yield B = 20.0^2:7 12'.4 MeV, corresponding to the g 
and /s uncertainty respectively. 

It is also interesting to notice that Ref. [27| predicts a 
total of two 1++ states in the bottom sector: a shallow 
bound state just below the BB* threshold and a deeper 
one with a binding energy about 140 MeV. The reason 
for the appearance of the two bound states is a unusu- 
ally big value of the BB*tt coupling. Within the present 
framework, and taking g — 0.55, a second bound state 
does not appear up to Ab — 1.2 GeV: for this cut-off, 
there is a low-lying bound state at B ~ and a second 
bound state at i? ~ 400 MeV. If we take g = 0.7, which is 
a relatively large value for the BBn coupling, the appear- 
ance of the second bound state happens at Ab — 1 GeV 
with a binding energy B ~ 175 MeV, not far away from 
the results of Ref. 27]. From the aforementioned values 
the possibility of a second bound state in the '^S'l — '^Di 
channel seems quite unlikely, in view of the values of the 
cut-off obtained in the present work. However, the re- 
quired cut-off and the position of a second bound state 
may depend on higher order contributions to the poten- 
tial. 



2. The Negative C-parity ^Si - ^_Di Bound State 

In second place we consider the negative C-parity part- 
ner of the previous state, that is, a BB* bound state 
in the '^Si — ^Di channel with total isospin 7 = and 
1^ . The complete wave function reads in this 



case 



BB* 
{pjm] 



(fc) 



D*0\ 



IB^B 



\B+B* 



\B*"B") 



\B*+B- 



(109) 



(pjm) (^) : 

(110) 

where p = +1, j = I and {l}^ — 0, 2. The spatial wave 

function, 5'(pj„j) (fc), is normalized and decomposed into 
an S- and D-wave component in exactly the same way as 
its positive C-parity partner. The short range potential 
takes the form given by Eq. (|108p and is subjected to 
the saturation condition, see Eq. (p^ . with fs = 195 ± 
10 MeV. The OPE potential has a overall minus sign 
with respect to the C = +1 case. 
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If we employ a sharp cut-off regulator with = 
3961^^ MeV, we obtain a binding energy of i? = 



6-4^46 MeV, where the error is again dominated by the 
l/mg corrections (i.e. the cut-off uncertainty) . However, 
contrary to the positive C-parity state, we do not have an 
external check of the validity of the saturation condition 
employed in this channel. This represent an additional 
source of error which cannot be easily estimated and we 
should therefore expect a much larger uncertainty in the 
position of this state. 



3. The ^Po Resonant State 

The last state we will consider in detail is a isoscalar 
(7 = 0), positive C-parity, ^Pq resonant state with J^'~^ — 
^ . For this state the wave function can be written as 



BB* 
{pjm) 



(k) 



\B+B*-) 



\B*"B°) 
- \B*+B- 



(pjm) 

(111) 



where p = —1, {l}^ = 1, j = and to = 0. Even though 
the intrinsic C-parity of this state is negative, the total 
C-parity is positive, as interchange of the meson and the 
antimeson will generate an additional minus sign in the 
wave function. The wave function only contains a P- 
wave component, which means that the separation of the 
radial and angular piece is trivial 



(112) 



^'(-ioo)(fc) = V47r*p(fc) 3^100 (fc) 
Finally, the normalization condition is 

"|vl/p(fc)|2^1, 



27r2 



(113) 



which only applies when the state is bound. 

For obtaining the resonant (bound) state energy we 
solve a one channel Lippmann-Schwinger equation with 
the OPE potential in the second (first) Riemann sheet 
(see Appendix [K\ for details). The potential is regulated 
with a sharp cut-off function, where the cut-off window 
is the one determined in the previous section (i.e. Ab = 
396t^|MeV). With this cut-off and g = 0.55 ± 0.1, we 
obtain a resonant state energy of 



Er 



0.7^ 



i.r 



■MeV, 



(114) 



where the upper error is to be interpreted as the disap- 
pearance of the resonant state, which eventually happens 
for the smaller values of the cut-off. On the other hand, 
the lower bound is compatible with a shallow bound 
state. The interesting feature of the theoretical ^Pq state 
is that, at leading order, it does not depend on the value 
of a contact operator, only on the estimation for the natu- 
ral value of the cut-off. Although the resonance vanishes 
for small values of g or of the cut-off, the existence of 



this states looks more probable than not. In addition, 
owing to the intrinsic negative C-parity of this state, the 
wave function can contain an appreciable B*B* compo- 
nent: this effect will reduce the energy of the system by 
a small amount, helping the formation of a bound state 
in this channel. 



4- The Isovector States and the Zi,(10610) Resonance 

The recent discovery by the Belle collaboration of two 
new resonances, the ^6(10610) and Zb(10650) [H, which 
lie a few MeV above the BB* and B*B* thresholds re- 
spectively, provides two interesting candidates for heavy 
meson molecular states. The quantum numbers of the 
two Zb states are I'^{J^) — 1"'"(1"'"). In the particu- 
lar case of the ^{,(10610), it is natural to interpret this 
resonance as a low-lying S-wave BB* state. Several the- 
oretical works 6_4-^] have appeared recently trying to 
explain the nature and properties of these two states. 

The S-wave molecular interpretation is however not 
completely trivial owing to the location of the Zfc(10610) 
at 4 ± 2 MeV above the BB* threshold. Within a po- 
tential description this requires an explanation in terms 
of a BB* resonance. In this regard, the existence of a 
resonant state in non-relativistic scattering depends on 
two ingredients, namely (i) a repulsive potential barrier 
at long distances and (ii) some sort of short range attrac- 
tion. If we consider the form of the OPE potential in the 
I — \, C = —1, ^Si — ^Di channel, we can appreciate 
that the OPE potential between two S-wave mesons is 
weakly repulsive and given by 



V;s(j=i)(r) = Tyc(r). 



(115) 



see Eq. ([89)1 . The shorter range attraction is provided by 
the S- and D-wave mixing, which in fact is able to bind 
the system for large enough values of the cut-off, see Ta- 
ble [ni In addition, being an intrinsic negative C-parity 
state, the BB* and B*B* channels can mix. If we take 
into account that the potential barrier provided by OPE 
is only about 1 MeV high at r = 0.5 fm, which is insuf- 
ficient for reaching the current position of the Zfc(10610) 
resonance, it is apparent that a molecular interpretation 
of this resonant state requires a coupled channel approach 
in which the B*B* is explicitly included. In principle the 
argument above certainly renders difficult the direct ex- 
trapolation of the (purely) contact theories employed to 
describe the Ar(3872) to the newly discovered Zf,(10610), 
as pion exchanges and coupled channels are expected 
to be important. However, the theoretical analysis of 
T(5S') -> /ibTT+TT- and hb{2P)Tr+n- decays_of Ref. [7l| 
suggests that the position of the eventual BB* state cor- 
responding to the Zh does not need to coincide with the 
position obtained from the Breit-Wigner parametrization 
employed in Ref. [l3|. In particular the .^(,(10610) can 
be located below threshold, around 5 MeV according to 
Ref. [7l|, reopening the possibility of a bound state in- 
terpretation. 
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number This will generate a contribution to the 
contact interaction in the isovector channel, Co = —6/fg, 
thus providing the missing attraction needed to generate 
a resonant (or bound) state in the ^(,(10610) channel. 



TABLE IV. Dependence of the binding energy of the two 
BB* isoscalar bound states with respect to different regu- 
lator choices. In particular we consider gaussian regulators 
with n = 1,2,3,4, labeled i?!"-^'^'^'*' and the sharp cut-off 
regulator B^c- As can be appreciated the central value of the 
binding energy can change moderately from one regulator to 
another. However, these variations are smaller than the un- 
certainty coming from other sources (in particular the 15% 
error associated with l/mg). 

In general, isovector heavy meson states are not usually 
predicted (or even considered) in potential models, the 
reason being that the expected strength of the potential 
is three times weaker in isovector than in isoscalar states. 
From Tablc|lT]we can see that the I{J^^) = 1(1"*"") state 
binds at a cut-off A ~ 1 GeV, which is definitively bigger 
than the obtained EFT cut-off. As ~ 0.4 GeV, or the 
values of the cut-off which bind the two isoscalar states, 
A ~ 0.4 and 0.7 GeV respectively E3. This does not mean 
however that the current estimation of the size of the cut- 
off is incorrect. A more natural explanation lies in the 
modification of the short-range dynamics. In particular, 
the saturation condition employed in the present work 
requires the contact interaction to be Co — —1/ f% in 
isoscalar channels and Co = in the isovectors, which 
penalizes the formation of / = 1 states. On the con- 
trary, if we assume a zero energy bound state in the 
Zb(10610) channel, the isovector counterterm needs to 
take the value Co = (—0.3 ±0.1)//^, which although 
non-zero, is still relatively small in comparison to the 
isoscalar counterterm. Moreover the size of the countert- 
erm is expected to decrease further if wc take into account 
that (i) the 2';,(10610) is above threshold (see however the 
previous discussion) and (ii) the mixing with the B*B* 
channel will provide additional attraction. In this regard, 
the only ingredient to accommodate the Zb(10610) state 
within the present framework is a small correction to the 
short range dynamics of the system in the line of 

(fc>cifc')=c^o^^*(|+^ 1;^) , (116) 

instead of the form given in Eq. (I38p , with 5 some small 



In a related note, the recent exploration of meson exchange 
in the BB* and B* B* systems of Ref. [TOj also indicates that 
the value of the monopolar cut-off which binds the Z(, (10610) 
channel is twice the size of the required cut-off for binding the 
I{J^^) = 0(1++) state, in qualitative agreement with the results 
of Table In] 



5. Higher Order Corrections 

In the previous calculations of the energies of the BB* 
bound states we have taken into account three sources of 
error: (i) the BB*tt coupling constant g, (ii) the weak de- 
cay constant fs and (iii) the uncertainty in the size of the 
1 /mg corrections for determining a suitable cut-off win- 
dow. An additional error source in the EFT formulation 
is the size of the contributions from higher order terms. 
From the EFT viewpoint, if an observable is computed 
at order P'^, the relative error for this observable is ex- 
pected to be 0{P''+^/Ao+^), where P and Ao are generic 
notation for the soft and high scales of the system. For 
a LO calculation (i.e. P°), the previous estimation will 
yield an expected error of ©(P/Aq). However, the EFT 
with pions and heavy mesons fields does not contain any 
correction at order P^ , the reason being parity conserva- 
tion. Therefore we expect the calculations of the binding 
energy to be accurate up to 

^EFT^^W^^l^^^ ^ (117) 

where B^^"^ is the full value of the binding energy in- 
cluding all the EFT corrections and P^^-* is the LO ap- 
proximation. In this context, P can be interpreted as 
the wave number of the bound state, that is P ~ 7^°^ = 
\/2 iipp, P(°), with iipp, the reduced mass of the heavy 
meson-antimeson system. This means in particular that 
the relative error is proportional to the binding energy; 
therefore, we can write 

PEFT ^ ^(0) ^ /^_^\ (118) 

where B^^^ is the maximum bound state energy which 
can be described by the EFT, which corresponds to 

Assuming a breakdown scale in the range Aq = 0.5 — 
1.0 GeV, we can estimate PJ'^'^ = 120 - 500 MeV for 
the DD* system and Pg^^'^ = 45 - 190M_eV for PP_*. 
As we are ignoring the mixing with the D*D* and B*B* 
channels, we expect the true breakdown scale to lie closer 
to Ao 0.5 GeV than 1 GeV. 



The model of Gamermann and Oset [2(| already contains this 
kind of corrections. However, their size may be far too small in 
the case of the BB* system, giving S ~ (mp/m~[)^ ~ 0.006. 
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As we can see, the lower estimation for the breakdown 
of the EFT description of the BB* bound states sug- 
gests moderate corrections to the estimation of the bind- 
ing energy of the I{J^^) — 0(1"'""'") state. This situation 
requires further analysis in order to check the reliability 
of of the results. In a standard cut-off EFT formulation, 
a practical way to estimate the size of the higher or- 
der contributions is to vary the cut-off within a sensible 
range E3- The underlying idea is that cut-off uncertain- 
ties are a higher order effect. In this sense, varying the 
cut-off mimics the effect of including (or excluding) the 
higher order contributions. However, in the non-standard 
formulation employed in this work, cut-off variations are 
used for estimating the size of l/mg corrections. In prin- 
ciple, this may be interpreted as the necessity of going 
to subleading orders to explicitly check the size of the 
higher order corrections. Nevertheless there is a second 
way of doing things which is to consider how the results 
vary with different regulators. In particular we can check 
the effect of using gaussian cut-offs of the type 

/(^)-e-^, (120) 

for different values of rt. As in the sharp cut-off case, we 
determine A by fixing the location of the AT (3872) state. 

The results for changing the regulator are shown in 
Table IIVI where we have considered the cases n = 
1,2,3,4[3- We have only considered the two S-wave 
isoscalar bound states, as the calculation of the P-wave 
resonant state energy is much more involved. As can 
be seen, the biggest change happens when the n — 1 
gaussian regulator is used, for which the central value 
of the cut-off is raised to A = 620 MeV and which gen- 
erates a change of 10 MeV in the binding energy of the 
0(1"^"^) ^Si - ^Di state (and 2.5 MeV in its negative C- 
parity partner). The other regulators generate however a 
much smaller change in the results. If we take the n = 1 
gaussian regulator as an upper bound of the LO uncer- 
tainties, we can appreciate that the EFT error is a bit 
smaller than the other error sources (in particular, the 
1/toq corrections). In addition, the range in which we 
expect the binding energy to lie does not change so much. 
The uncertainties follow the expectations of Eq. (jll8l) . 
that is, the relative error grows with the binding energy. 
On the other hand, the 40 — 50% relative error in the 
0(1"'""'") state is consistent with the lower estimations of 
the breakdown scale, indicating that we are well within 
the range of validity of the EFT with pions. 



See for example Ref. |41|| for the application of this idea in the 
context of nucleon-nucleon scattering 
^® It should be noticed that 2n must be higher than the order 
P'^ at which the computations are done to avoid contamination 
of the auxiliary cut-off scale at orders below that of the EFT 
calculation. For the LO calculation of the present work, this 
condition does not translate into a constraint for the value of n. 
However, if we go to order P^, we should employ at least n > 2. 



6. Further Uncertainties in the BB*n Coupling 

As previously mentioned in this section, a particular 
problem we encounter when considering the OPE poten- 
tial in the bottom sector is the determination of a suit- 
able value of the BB*n coupling. The value we have 
employed, g = 0.55 ± 0.10, approximately encompasses 
most results from lattice QCD, which usually range from 
g = 0.44 [Hi] to g = 0.63 [Hi 0. However, this value of 
the BB*Tr coupling does not reflect all the theoretical un- 
certainties involved in the determination of this quantity 
(see, for example, Refs. [5^ [tI] for a compilation of val- 
ues). In particular, smaller values of this coupling may 
be possible, even in the range g = 0.3 — 0.4, a case which 
we will consider here. A small BB*tt coupling translates 
into a weaker OPE potential, as its strength scales as g^, 
and a much weaker chiral TPE (oc g^), meaning that the 
higher order corrections owing to pion exchanges will be 
strongly suppressed. 

If we take g = O.ST^p gg as a reference value, which 
was obtained in Ref. [59j from the Dyson-Schwinger equa- 
tions of QCD, we predict the binding energies -6(1"'"'") = 
16l9^MeV and B{1+-) = 9.2^^;^ MeV for the positive 
and negative C-parity ^Si — '^Di isoscalar states respec- 
tively. These new values mostly overlap with the predic- 
tions corresponding tO(7 = 0.55±0.10, which is consis- 
tent with the observation that contact operators are the 
dominating mechanism in the formation of S-waves BB* 
bound states. On the contrary, the '^Pq (0 '") resonant 
state disappears, as it depends crucially on the strength 
of the OPE potential. For the isovector ^{,(10610) state, 
the results do no significantly change: the contact term 
still needs to be of the order of Co — —0.3//^ to bind 
the state. 



V. CONCLUSIONS 

We have considered the DD* and BB* two meson sys- 
tems within the framework of a pionless and a pionfuU 
(or chiral) cut-off EFT at LO. In the charmed sector, the 
existence of the A(3872) state, together with the satura- 
tion hypothesis for the low energy constant Cq^ , sets 
the conditions for the applicability of the EFT formula- 
tion employed in this work and results in a natural value 
of the cut-off Ax ■ The determination of the low energy 
constant and the value of the cut-off in the BB* case 
requires invoking HQS to overcome the absence of exper- 
imental information for this two body system. In this 
regard HQS is able to correlate the charm and bottom 
sectors, and, in addition, provides error estimates for the 



It should be noted that in most of lattice QCD calculations the 
value of the BB*it coupling is computed in the mq — >■ oo limit, 
probably with the exception of Ref. 55], in which g = 0.58 is 
obtained. In general, we expect g{mQ = nii,) > g{mQ — > oo). 
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contact term and the cut-ofF in the BB* system. This 
in turn allows to assign errors to the resulting binding 
energies of the possible BB* bound states. 

The present framework predicts the existence of three 
isoscalar BB* states with positive C-parity: two — 
'^Di states with positive and negative C-parity with a 
binding energy of 20 MeV and 6 MeV respectively, and a 
'^Pq (0 ^) resonant state which lies almost at the BB* 
threshold. The different error sources result in relatively 
large uncertainties in the previous estimations. In ad- 
dition, we expect moderate corrections from subleading 
order contributions to the chiral potential between the 
two heavy mesons. However, the higher order correc- 
tions are probably smaller than the current uncertainties 
stemming from the use of the approximate HQS. Never- 
theless, the existence of the bound states is a conclusion 
that will likely remain unchanged: if the binding energy 
is lowered, the subleading order corrections are expected 
to decrease, resulting in a stabilization of the results. 

The applicability of the present approach to negative 
C-parity states is not entirely free of problems, one of 
them being that there is no direct experimental evidence 
supporting the short-range dynamics employed in the 
present work. This feature translates into an additional 
error source in the isocalar C = — 1 "^Si — '^Di bound 
state which we are not able to estimate. In addition, the 
existence of the Zb(10610) state points out to small but 
significant deviations from the form of the short range 
interaction employed in the present work. A second is- 
sue with the intrinsic negative C-parity states is that they 
tend to couple more strongly with the B*B* system than 
their positive C-parity counterparts, an effect that is ex- 
pected to increase the overall attraction of the system 
and help the formation of bound states, specifically the 
aforementioned Z(,(10610). 

The extension of the present framework to higher or- 
ders in the EFT formulation, or to other heavy meson 
systems, such as PP and p*p* ^ is left for future research. 
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Appendix A: Resonant State Equation 

In this section we extend the Lippmann-Schwinger 
equation to the second Riemann sheet for finding res- 
onant and virtual states. For simplicity we will work in 
the partial wave decomposition corresponding to the cen- 
tral potential case; the extension to coupled channels and 



tensor forces is straightforward. The starting point is to 
use the vertex function, instead of the more usual wave 
function. The vertex function is related with the residue 
of the T-matrix at the pole energy, that is 

hin [E - EB){k\Ti{E)\k') = ^B,i{k) (l>B.i{k') , 

E^Eb 

(Al) 

consequently the relation between the vertex and the 
wave function is given by 

= Go0s(fc). (A2) 

Inserting this relationship into the bound state equation, 
Eq. (fT3)) . and extending the equation to arbitrary ener- 
gies, we obtain 

(t)B,iik) = 

--^J ^i^miq)mBAk'), (A3) 

where = —7^ = 2/ii5. The advantage of the vertex 
equation is that it contains the resolvent operator within 
the integral, which makes it possible to select the Rie- 
mann sheet by deforming the integration contour around 
the k^ = singularity, see Ref. [75;] . The vertex equation 
only has solutions for the energies at which the T-matrix 
has a pole. 

In principle, to find the position of resonant and virtual 
states it is enough to analytically extend the solutions of 
the vertex equation to the second Riemann sheet. How- 
ever, in numerical calculations the previous is not triv- 
ial. The numerical evaluation of the resolvent operator, 
Go{E), always chooses the first Riemann sheet, prevent- 
ing us from finding either resonant or virtual states. The 
solution is to force the selection of the second Riemann 
sheet, for example, by changing the integration contour 
as previously commented. Here we will use instead a 
more informal derivation for extending Eq. (jA3p to the 
second Riemann sheet. We will assume that the physi- 
cal scattering region corresponds to E + ie, with E real 
and positive. If we consider E — ie instead, the resolvent 
operator takes the form 

-^=Vi-^)+^7^^S{k-q), (A4) 

where V denotes that the principal value should be taken 

2 

and with E = However, if we were in the second 

2/i ... 
Riemann sheet, we would need the imaginary piece of the 

resolvent operator to be negative instead of positive: we 

are moving from E -\- ie to E — ie in a, continuous manner, 

which means that the imaginary piece should be the same 

in i? -f ie as in £' — ie. A practical solution is to add the 

imaginary piece directly into the resolvent operator, that 

is, in the second Riemann sheet we substitute the original 

resolvent operator Gq^"* by a new resolvent operator Gq^^^ 

gII\e) ^ Gi"^ = Gi^^E) - i2TT ^ 5{k - q) . (A5) 
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From the point of view of the bound state equation, the 
previous changes amounts to the substitution 



1 



k'^ - ff2 fc'2 _ q2 q 



^ +i-6{k'-q), (A6) 



within the integral in Eq. (jASp . This change leads to the 
following set of equations 



^ f k'^dk' 



n^J k'^-q^ 



(fc|y,|fc')0i'.P(fc'), (A7) 



^ f k'^dk' 



n^J k'^-q^ 



{qmk'}<j>^^^,\k'), (A8) 



which are equivalent to the analogous set of equations 
obtained in Ref. [t^. Rearranging the different terms in 

the previous equation to eliminate the cj)^j^\^ {q) term, we 
arrive at 



^ / k''^dk' ^^\nj\A")fl '\ 



k'^ - q' 



where Wi{q) is defined as 



■{k\WM)W)<I^Vj{k'), (A9) 



{k\Wi{q)\k') = {k\Vi\k')-l 



fiq {k\Vi\q) {q\Vi\k') 
TT l + ^if{q\V^\q) 



(AlO) 



Depending on the value of q'^, that is, of E, the pre- 
vious equation will look for resonances (Re(i?) > and 
Im(£') < 0) or virtual states (Re(£') < and lm{E) ^ 0). 
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